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Accuracy of a relativistic weak-coupling expansion procedure for solving the Hamiltonian bound- 
state eigenvalue problem in theories with asymptotic freedom is measured using a well-known matrix 
model. The model is exactly soluble and simple enough to study the method up to sixth order in 
the expansion. The procedure is found in this case to match the precision of the best available 
benchmark method of the altered Wegner flow equation, reaching the accuracy of a few percent. 
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I. INTRODUCTION 

This article describes a foretaste study of accuracy for 
a recently proposed Hamiltonian weak-coupling expan- 
sion procedure that in principle can start from a local 
asymptotically-free quantum field theory and produce 
sufficiently small relativistic effective Hamiltonian eigen- 
value problems that may be soluble on a computer and 
yield the wave functions of bound states in that the- 
ory. The study is needed to determine the prospects of 
reaching a reasonable accuracy in the expansion since the 
strong coupling constant rises when the renormalization 
group scale is lowered toward the scale of binding mech- 
anism 0,13 ■ So far, successful approaches required dif- 
ferent formulations of the theory at the bound-state and 
high-energy scales, such as Wilson's lattice and Feyn- 
man's diagrammatic techniques H Q], Hi E| . They also 
used approximations such as the non-relativistic limit 
that applies in the case of heavy quarkonia 0, IE EJ 
and helps in constructing effective theories in analogy 
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The approach discussed here is the renormalization 
group procedure for effective particles (RGEP) [T^ that 
stems from the application of the similarity renormaliza- 
tion group procedure |l5j| to the light-front Hamiltonian 
of QCD In one and the same scheme, RGEP pro- 

duces the asymptotically free coupling constant in the 
Hamiltonians for quarks and gluons [18j, provides the 
conceptual framework for constructing the whole renor- 
malized Poincare algebra in terms of the creation and an- 
nihilation operators for effective particles 0, and leads 
to a simple first approximation in the case of heavy 
quarkonia poj . The procedure is boost invariant and 
raises a hope for connecting the constituent model for 
hadrons at rest plj and the parton picture in the infi- 
nite momentum frame The key ingredient of the 
procedure is the vertex form factor / that multiplies all 
interaction terms. It falls off quickly to zero when the 
change of an invariant mass in an interaction vertex ex- 
ceeds the width parameter A. This width variable is also 
the renormalization-group evolution parameter. The ef- 
fective particles that correspond to a small value of A can- 
not be copiously created because the form factor makes 



the interactions effectively weak even if the coupling con- 
stant becomes large. The internal structure of the effec- 
tive particles at such small A may still be given by the 
parton-like picture in terms of constituents that corre- 
spond to a much larger value of A ~ Q in the same flow, 
where Q is the momentum scale of the external probe 

M 

Qualitative accuracy studies of a similarity scheme 
were performed before [24[ using the elegant Wegner flow 
equation [2 a, |26( , which was invented independently for 
solving Hamiltonian eigenvalue problems in solid state 
physics. Unfortunately, it was found that the Wegner 
equation was not suitable for a straightforward weak- 
coupling expansion beyond second order. The coeffi- 
cients of the expansion grew too fast and alternated in 
sign, which lead to erratic results for bound state ener- 
gies with no signs of convergence. But the Wegner equa- 
tion can be modified within the similarity scheme |27j 
and the improved equation provides the benchmark here 
for estimating the accuracy prospects for the relativis- 
tic RGEP procedure. The comparison with the altered 
Wegner equation and the fact that improvements are also 
needed in condensed matter physics [2^, HtJ |3(j , imply 
that the effective particle approach may also find applica- 
tion outside QCD, i.e., wherever the dynamical couplings 
increase in the flow of Hamiltonians toward the region of 
physical interest. 

The RGEP strategy that is tested here is to start with 
a regulated H of the theory to be solved (this H provides 
an initial condition for the differential equations of RGEP 
at A = oo), and to evaluate H\ with A on the order of 
a bound-state energy as a series in powers of a coupling 
constant g\. After evaluation of H\, one calculates its 
matrix elements in the basis defined by eigenstates of the 
Hamiltonian Hq\, which is obtaind from H\ by setting 
gx = 0; H \\n) = E n \n) for all values of the label n. 
Suppose that the labels are ordered so that E m < E n 
implies m < n and the initial Hamiltonian is regulated by 
forcing matrix elements H mn = (m\H\n) to vanish unless 
M < rn,n < N with certain ultraviolet cutoff number N 
and infrared cutoff M . The RGEP procedure is designed 
in such a way that the matrix elements (m\H\\n) quickly 
tend to when \E m — E n \ grows above A, see Fig. ^ 
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FIG. 1: This out-of-scale figure illustrates Eq. Q for the 
matrix W\ mn with indices M < m,n < N . The matrix of 
H\ is represented by the large square. Matrix elements of H\ 
outside the diagonal band of width A are equivalent to zero. 
In the model study: E N ~ 65000 GeV, E$ ~ 4 GeV, and the 
width A ~ 2 GeV. 

The next step is to focus attention on the window W\ of 
matrix elements of H\ among the basis states that have 
energies similar to the energy scale E of the physical 
problem at hand, i.e., 

W X mn = {m\H x \n), (1) 

with M < m,n < N and E Al < E < E N . Since only 
states within the width A on the energy scale can di- 
rectly interact with each other, states that differ in energy 
from E by much more than A are usually not important 
pil W\ . They can be important as long as the coupling 
strength can overcome the difference in energy, but it 
does not matter here because the coupling constant is 
assumed to not grow to large values. The next step is 
to solve the non-perturbative eigenvalue equation for the 
matrix W\ mn by diagonalizing it on a computer. The 
middle-size eigenvalues of W are expected to be close to 
the exact solutions with accuracy that depends on many 
factors in the procedure. These dependencies need to be 
estimated. The main question addressed here concerns 
the accuracy of the weak-coupling expansion for W\. 

An exact RGEP procedure provides W\ mn whose 
structure depends on A but the eigenvalues of the middle 
size in the window spectrum do not. Becasue the cou- 
plings to the states outside the window are ignored, the 
eigenvalues of sizes at the limits of the window spectrum 
cannot be accurate even if the window is calculated ex- 
actly. Now, when one calculates W\ in the expansion 
in powers of g\ to some order, and extrapolates the re- 
sult to g\ ~ 1, considerable errors can ensue because 
of the missing terms in the series. Moreover, one never 
knows the right value of g\ at given A from the theory. 
Therefore, one must fit g\ to the bound-state observables 
and perform consistency checks for a whole set of them. 



The smaller is A the larger is g\ and the more signifi- 
cant are the errors of perturbation theory in evaluation 
of W\. But at the same time the larger is A, the larger 
must be the range [M, N] of basis states needed in the 
computer diagonalization of W\ . A compromise must be 
made and the critical question is how close one can come 
to a true solution using RGEP equations. The question 
is essential for the prospects of applying RGEP to QCD. 
A well-known asymptotically free matrix model is used 
here to find out what level of accuracy can in principle 
be expected. One should stress at this point that even 
if the test gave a promising result in the model, the util- 
ity of the procedure would remain only tentative until 
the actual calculations in realistic theories are performed 
and display signs of stability as functions of the order of 
the expansion, size of A, variation of the window, and 
other more specific features of the eigenvalue problems 
at hand. 

Section [D] defines the test model and quotes equations 
used for calculating H\ in the benchmark of the altered 
Wegner 's flow and in the RGEP case. Section [Till de- 
scribes results obtained using six different ways of fitting 
the coupling constant to the exact spectrum. These ways 
are labeled throughout the paper by letters A, B, C, D, 
E, and F. Each of these versions is studied in six suc- 
cessive orders of the weak-coupling expansion. Section 
IIVI provides a brief conclusion. Appendix lAl contains the 
generic analytic RGEP formula for 4th order calculations 
that applies to arbitrary iJ, and Appendix ]B\ provides 
the numbers that illustrate in detail what happens in the 
expansions including from one to six orders. 



II. MODEL 

The matrix elements of the model Hamiltonian used 
here for estimating the accuracy of RGEP are [24L [27j 



Elmn — E n 8 mn g \l E ra E n , (2) 

where E n = b n , b > 1, and n is an integer, with a con- 
vention that the energy equal 1 corresponds to 1 GeV. 
The model diverges and needs to be regulated by the 
ultraviolet cutoff A = b N that limits the allowed ener- 
gies E n from above. Also a low energy cutoff, b M , is 
introduced for making the Hamiltonian matrix finite and 
enable exact computations, but the lower bound is of no 
physical consequence for the results reported here. Thus, 
the subscripts in Eq. J2J are limited to the range [M, N] , 
M being large negative and N large positive. The ul- 
traviolet renormalizability of the model, its asymptotic 
freedom, and its lack of sensitivity to the infrared cutoff, 
were described in |24L |27| . 

Two values of the cutoff N and the corresponding cou- 
pling constant g — gN are used in this study of RGEP: 
gie = 0.060600631 and g 20 = 0.04878048667. The two 
cutoffs are introduced to verify the accuracy of renormal- 
izability in the RGEP scheme. It is known that for such 
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large values of N the effective dynamics with A ~ 1 is 
practically independent of N in the case of altered Weg- 
ner's equation, and one can verify how well the RGEP 
approach satisfies this condition. At the same time, this 
condition puts constraints on the range of changes that 
one can consider in varying the RGEP generator without 
interference with renormalizability (see below) . The cou- 
pling constants gi$ amd 520 are fitted to obtain the exact 
bound-state eigenvalue E = — 1 with 8 digits of accuracy 
for b = 2 and M = —21, as in 

El- 

With these choices, 
the Hamiltonian H is a (N - M + 1) x (N — M + l) ma- 
trix with N — M positive eigenvalues, and one negative, 
spanning the range of energies between Em ~ 0.5 KeV 
and E N ~ 65 or 1000 TeV. 

The similarity renormalization group procedure for 
Hamiltonians that leads to the altered Wegner equation 
pBll^ and provides the benchmark here, can be written 
in the differcntail form as 



dU 
~dX 



[F{H},H] 



(3) 



with the initial condition TL = H when A = 00. The 
initial condition contains counterterms but the similarity 
analysis showed 0] that the structure of the countert- 
erms is simple and the presence of them is equivalent for 
large N to making g in H depend on N. This is precisely 
what is done by the fitting mentioned earlier that guar- 
antees that the eigenvalue E = — 1 stays unchanged for 
different TVs. All small eigenvalues are then independent 
of N. The generator of the similarity transformation can 
be written as (T> m = Ti. mm ) 



(m\F{H}\n) = h mn (V m - V n )H r , 



(4) 



Different choices of h mn lead to different matrix elements 
of the renormalized Hamiltonians. Assuming 

, _ , ds_ 

"run — fymn 7 -. j 

dX 

one obtains Wegner's equation when <p mn = 1, and s = 
1/A 2 plays the role of the original Wegner parameter I 
|25| . The altered Wegner equation has [27j 



(5) 



1 



1 + \m — n\ ' 



(6) 



and this new equation is referred to as the benchmark. 

In the plain perturbative RGEP procedure, the matrix 
elements of H\ in the effective basis states associated 
with the scale A, are obtained from the matrix elements 
of an auxiliary Hamiltonian Ti in the initial basis |l7j| . 
The structure of Ti. is given by 



Tt = TLq + fGi ■ 



(7) 



The Hamiltonian 7i is equal to the initial H , which is 
the free part of H, i.e., H with g = 0. / denotes the form 
factor that can be written in all matrix elements using 
eigenvalues of Hq , and reads 



fnui — exp 



(En 



E„ 



A 2 



(8) 



The RGEP equation for Qi is 



dGj_ 
dX 



(9) 



where the curly bracket around an operator has the fol- 
lowing meaning in terms of the matrix elements, 



E n — E n 



(10) 



A priori, the optimal choice for cf> mn in the RGEP could 
be different from the one that optimized the benchmark 
|27) . But we have studied various factors (j) mn a la [53 
and found that c ~ 1 is also the best choice to make 
in RGEP, for similar reasons. In addition, it is useful 
for the test that these factors are made equal since then 
the first order calculations give identical results in both 
approaches. The factor <j) mn is included in the analytic 
4th order RGEP formulae provided in Appendix IA1 

The RGEP Eq. 10 cannot be integrated exactly on 
a computer as easily as Wegner's equation can, because 
it contains the derivative of Qi on its right-hand side. 
One has to solve a complex linear problem to extract Q\. 
But the altered Wegner equation provides a perturbative 
benchmark pattern that is known to approximate an ex- 
act solution well and one can estimate the accuracy of 
RGEP by comparison. 

In the perturbative evaluation of H\, the RGEP cal- 
culus is free from the problem of extracting Q\ because 
the derivative is computed order by order and all terms 
needed on the right-hand side are known at each succes- 
sive order from the lower order results. In fact, the new 
procedure is designed for a perturbative approach. It is 
simpler than the Wegner case in the sense that the form 
factors guarantee the band structure in perturbation the- 
ory so that this structure does not have to be recovered 
from and controlled in the evolution of specific matrix 
elements. Small energy denominators that might other- 
wise lead to infrared singularities are excluded by design 
of RGEP. It also does not generate any terms with in- 
verse powers of A in the coefficients of products of the 
interaction Hamiltonian (see Appendix 1X1 and Eq. (3.2) 
in |27jp. Such terms can in principle lead to a variety 
of mixing effects in the evolution of H\, when one uses 
the Wegner equation. On top of these purely quantum 
mechanical features, the perturbative RGEP calculus is 
capable of respecting seven Poincare symmetries in an 
economic way and has a potential to obtain the remain- 
ing three [ljj- At the same time, it preserves the cluster 
decomposition property \?>\\ in the effective interactions. 
All these features are desired for the description of rela- 
tivistic particles using theories with asymptotic freedom. 
But the accuracy of the weak-coupling expansion for win- 
dows W\ in RGEP must be measured against the bench- 
mark to estimate the cost of the apparent advantages in 
terms of precision. 
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III. ACCURACY OF RGEP 

The results of weak-coupling expansions up to the first- 
order terms in the RGEP and altered Wegner equation 
are identical, and read 



approximates the exact one as closely as possible. We 
display results for six options that illustrate the depen- 
dence of results on such choices, labeled by A, B, C, D, 
E, and F. All methods are based on the minimalization 
of certain function K(g\). We use 



H 



g\ V E m E n fn 



(ii) 



The form factor f mn causes that the interaction Hamilto- 
nian matrix is narrow on the energy scale and has width 
A, see Figure ^ It is clear that in the case of Eq. 111(1 
the effective coupling constant can be extracted from the 
matrix TL using the formula 



9\ 



1 



Wm,m(A) 
Em 



(12) 



which is analogous to the Thomson limit in QED when 
M is large and negative. The same Eq. lfrl|) is used for 
defining g\ also in higher order calculations. 

In the weak-coupling expansion in powers of g\ to order 
k (see e.g. 0), 



H = Ho + g\Hi 



g\n2 



(13) 



Appendix ^ contains analytic expressions for Jik with 
k = 1, 2, 3, and 4. The right value of g\ can be found 
by solving the flow equations with the initial condition 
9oc = gN exactly. Then, one can check the accuracy of a 
perturbatively computed H\ by diagonalizing it for the 
exact value of g\, and by comparing the resulting bound- 
state eigenvalue with E = — 1 that was secured to exist 
by the initial choice of g^ The accuracy test for 

RGEP is carried out differently because the exact solu- 
tion of Eq. © is not known. This situation is analogous 
to QCD where one can use perturbation theory to cal- 
culate H\ but an exact value of g\ is not available as 
a function of gN- For the purpose of the accuracy test, 
the exact spectrum of the model is treated as experimen- 
tal data. An approximate value of g\ is found by fitting 
some eigenvalue of the perturbatively calculated W\ to 
the data, or by performing a fit for a whole group of 
eigenvalues. Then one checks how well the bound-state 
eigenvalue E = — 1 is reproduced using the best-fit value 
for g\. In principle, one could fit g\ at one value of 
A that is most convenient for that purpose, evolve this 
value using RGEP to the new A that is most suitable for 
the bound-state calculation, and then compare the cal- 
culated spectrum with data [E H2i US ■ l n f act ' tne 
model used here can be used for testing accuracy of such 
procedures in a comprehensive way. This type of tests 
may help in narrowing the current spread of estimates 
for a s that come from various sources pi) , by distin- 
guishing theoretical procedures of least ambiguity. But 
the accuracy of RGEP is checked here using one and the 
same scale for fitting g\ and calculating the bound-state 
energy, for simplicity. The scale we choose is A = 2. 

There exist infinitely many options for how one can 
fit g\ so that the spectrum of the window W\ at A = 2 



and 



K s (g x ) = Zj2 



Vn ~ Vn+1 
n n+l 



- 1 



(14) 



(15) 



where v en is the exact eigenvalue of number n, and v n is 
the corresponding eigenvalue of W\ when W\ is derived 
in a given order k. Z is the normalization constant equal 
to the inverse of the number of terms in the sum. Z is 
not important in the minimalization of K(g\) as function 
of g\ . The subscript r refers to the ratios of eigenvalues 
used in K r , and the subscript s refers to the splittings 
between the eigenvalues used in K s . The six choices differ 
by which function K(g\) is used and what is the range 
of summation over n in Eqs. I|14|) or (|15fl . 

The eigenvalues of H are numbered from M to N in the 
order in which they appear on the diagonal of H\ when 
A — > in the benchmark calculation. This numbering is 
also applied to the corresponding rows and columns of 
the matrix H\ mn . The window W\ is always chosen to 
extend from M = -8 to N = 2 [23, and it has 11 eigen- 
values. These are numbered in the same order as for H\, 
with the bound-state eigenvalue having number n . We 
distinguish two eigenvalues that are closest in modulus to 
the bound-state energy E = — 1, one just smaller than 1, 
with number n s , and one just larger than 1 with number 
ni- The six fitting procedures are designated as follows: 



(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



In cases E and F, the two smallest and two largest eigen- 
values are dropped because they are too much distorted 
by the edge of W\, as explained in Section |I| The de- 
scription of results obtained in the benchmark and tested 
approach in six successive orders of perturbation theory 
in each of these fitting schemes involves 72 results. The 
benchmark results are labeled "Wegner" and the tested 
case is labeled "RGEP." Appendix iBl contains all per- 
tinent numbers. An example is given below to help in 
reading the figures and tables in the appendix. Other- 
wise, only main features of the results are explained. 

Fig. |2 corresponds to the case C of Eq. l(T%|) . The 
shape of K(g\) clearly selects the value of g\ preferred by 
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FIG. 2: Plots of K(g x ) from Eq. (0), in the altered Wegner 
case (benchmark) , and in the RGEP procedure, as marked, in 
the fit C of Eq. 118H and six successive orders of perturbation 
theory used in evaluation of the window W\, with the renor- 
malization group parameter A = 2 GeV, as functions of the 
a priori unknown value of the coupling constant g\ at this A. 
K(g\) measures deviation of two selected eigenvalues of the 
window W\ (those which are nearst-in-size to the bound-state 
eigenvalue) from their exact counterparts. The shapes of the 
functions point toward the required values of g\. This exam- 
ple is provided for explanation of how to read other examples 
given in Appendix IBl 



a given fitting procedure. Similar plots in other cases are 
given in Fig. IBJin Appendix iBl The numbers that result 
from the fits for the coupling constants, along with the 
corresponding bound-state eigenvalues of W\, are also 
tabulated in Appendix IBl 

The summary of results for the coupling constants ob- 
tained in the fits is given in Fig. [3] The numbers 1 to 6 on 
the horizontal axis correspond to the order of perturba- 
tion theory in the evaluation of W\. The columns labeled 
A to F correspond to the algorithms given in Eqs. Hltj|) 
to (|2ip. The RGEP calculation equals the benchmark in 
the first order results. It is also visible that the fits of g\ 
in the benchmark case consistently reproduce the exact 
value of g x = 0.2852 at A = 2 GeV H3- The RGEP 
displays similar stability and coalesces around 0.3. How- 
ever, when the nearst neighbor level with energy larger 
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FIG. 3: Coupling constants g\ obtained from fits described in 
Eqs. 1 Kit to (112 1 ft . gathered in the columns labeled by A to F, 
in the benchmark case (Wegner) and the tested case (RGEP). 
The six entries in each column correspond to six successive 
orders of perturbation theory used in the calculation of the 
window W\ with A = 2 GeV. 



than 1 is included in a fit as the only one (case B), or 
together with a nearest lower level using ratios of eigen- 
values (case C), or together with the nearest lower level 
but using ratios of splittings between the two eigenval- 
ues (case D), a visible variation in the fits occurs. This 
variation can be attributed to the lack of accuracy in the 
calculation of the nearst higher level, since in the cases E 
and F that include additional five lower levels, the higher 
level becoming much less significant, the fits resemble 
case A where the higher level is absent. This result sug- 
gests a rule for fits in future calculations that they should 
be focused on states with eigenvalues as far as possible 
from the bounds implied by the window choice and A, in 
order to avoid the lack of convergence. Note that even 
in the benchmark case the 4th order calculation has to 
be corrected in the orders 5th and 6th to bring the accu- 
racy into the few percent range around the exact value 
of g\. The same effect is observed in the bound-state 
eigenvalues themselves. 

The summary of results for the bound-state eigenval- 
ues is provided in Fig. 0]in a one-to-one correspondence 
to Fig. |3| The eigenvalues are obtained by diagonaliza- 
tion of windows W\ that are calculated in the six con- 
secutive orders of perturbation theory (indicated on the 
horizontal axis in the individual columns), using the cor- 
responding values of g\ from Fig. [3] (the numbers are 
given in Appendix lB)l . It is visible that the RGEP pro- 
cedure matches accuracy of the altered Wegner equation. 
One also sees the dramatic effect of the attempts to in- 
clude the next higher level. The results clearly point 
out that both procedures should not be used for calcula- 
tions of energy levels close to the size of A. The fourth 
order calculations achieve accuracy on the order of 5%. 
This is encouraging, because orders 5 and 6 lead to even 
better results and one can expect that fits for g\ can be 
improved by focusing on the properties of low energy lev- 



6 



Bound state energy 



■ - Wegner 
t - RGEP 



E 

-0.4 



A 


B 


C 


D 


E 


F 








I V 








■ 












\ ▼ 




i V 
» \ 










* 






















it 










■ ■ 

. ■ 




y 


if 




■ 










— ■ 




123456 


123456 


123456 


123456 


123456 


123456 



2nd order 



3rd order 



Order 

FIG. 4: Energies of bound states obtained from a non- 
perturbative diagonalization of window matrices W\ with 
A = 2 GeV. The matrices were derived in perturbation theory 
and evaluated using the corresponding values of the renormal- 
ized coupling constants, g\, from Fig. |3 The eigenvalues are 
displayed in the same convention and in the one-to-one cor- 
respondence to Fig. Numerical values are tabulated in 
Appendix iBl The exact result is -1. 



els, including properties other than just eigenvalues. It is 
hard to imagine that such focus could lead to worsening 
of the accuracy. 

The last question concerning accuracy of the RGEP 
method concerns renormalizability, which can be studied 
in analogy to |27| . However, one has to measure the sen- 
sitivity of the effective theory at the scale A to the cutoff 
A = b N in a whole range of coupling constants, because 
the exact value of g\ is not known in the RGEP proce- 
dure. Using giQ and 520 from Section^ we calculate W\ 
starting from H with the two different values of A for 
N = 16 and N = 20, and we compare the resulting ma- 
trix elements of W\ at A = 2 GeV. For 6 = 2 the cutoff 
is changed by the factor 16, i.e. from about 65 TeV to 
about 1000 TeV. The divergence in the bare theory is log- 
arithmic. The results could change at the rates implied 
by the change in the logarithm of A, i.e., about 25%, but 
in the presence of a proper set of counterterms |l5j | one 
expects no change to occur. Fig. [S] shows the measure 
of changes in W\ that are obtained here with only one 
counterterm, which amounts to the change of the cou- 
pling constant from gi§ to 520- The function plotted in 
ISJis defined as 
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(22) 



The results change with the order of perturbation the- 
ory used for evaluation of W\ . The sensitivity of results 
to the change of A in 1st order is the same in RGEP as 
in the benchmark case [27|. Inclusion of higher orders, 
from 2 to 6, exhibits slight variations to the advantage of 
one or the other method, as shown in Fig. (the solid 
line for the benchmark and the dashed one for RGEP). 
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FIG. 5: Plots of the renormalizability measure R(g\) from 
Eq. I22li for the RGEP procedure (dashed lines) and the al- 
tered Wegner equation (solid line) . The last diagram displays 
the ratios Rw / Rr (with Rw standing for the benchmark R, 
and Rr for the R in the RGEP procedure), and the integer 
labels indicate the corresponding orders of the weak-coupling 
expansion. 



The last diagram in Fig. [SJshows the corresponding ratio 
Rw/Rr, and it demonstrates that the renormalizability 
condition is satisfied with the same accuracy in both ap- 
proaches. A closer comparison requires that the bench- 
mark case (or Rw) is taken for g\ = 0.2852, and RGEP 
(or Rr) for some g\ ~ 0.3 ± 0.005 (see Appendix lB|l . but 
these corrections are negligible at the current stage of the 
analysis. 



IV. CONCLUSION 

The RGEP weak-coupling expansion achieves precision 
comparable with that of the best available benchmark 
method of the altered Wegner equation in the case of a 
simple matrix model for the Hamiltonian formulation of 
theories with asymptotic freedom and bound states. The 
few percent accuracy is reached by introducing a factor 
of <f> in Eq. ©, which is similar in both methods. Since 
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the RGEP procedure is designed for application to rela- 
tivistic quantum field theory, the model test calculation 
provides a comprehensive outline of steps that can be re- 
peated in realistic cases, especially in the theory of quarks 
and gluons. On the other hand, knowing that Wegner's 
equation is useful in condensed matter physics, one can 
expect that the altered Wegner equation and the RGEP 
procedure may also find applications in many-body the- 
ory. The model study indicates a possibility that the 
weak-coupling expansion may lead to a systematic ap- 
proximation scheme despite the growth of the coupling 
when the characteristic scale A of the effective theories is 
lowered. This is indicated in columns E and F in Figs. |3| 
and 0J where one can see the convergence of the coupling 
constant to a stable value and the corresponding con- 
vergence of the bound-state energy to the exact result. 
One should also note that already second order calcula- 
tions render effective window Hamiltonians W\ that can 
produce the bound-state eigenvalue with accuracy better 
than 10%. 



The fact that both methods have similar accuracy sug- 
gests that they already show the range of calculational 
power that is available through a plain expansion in a sin- 
gle coupling constant such as the one defined in Eq. (JT2J. 
In order to move beyond the 1% accuracy level one has 
to achieve better understanding of the structure of H\. 
For example, there exist terms in H\ with specific depen- 
dence on the eigenvalues of Hq and A, like (E m — E n ) 2 /X 2 
in the model. One may hope to understand how the coef- 
ficients of these terms depend on a more suitably defined 
coupling constant than the one given Eq. (|12|l . It is con- 
ceivable that such understanding may further improve 
the weak-coupling expansion. Note also that none of the 
renormalization group universality features were so far 
explicitly employed in the plain expansion tested here. 
Such options for improvement may depend on a theory. 
One has to study specific theories using concrete versions 
of the RGEP procedure in order to identify the dominant 
terms and their universal behavior. 



The accuracy study for the RGEP procedure shows 
also that one should be able to carry out similar tests 
for any other approach to the bound-state problem in 
asymptotically free theories. For such a test to become 
possible, the approach in question would have to be un- 
derstood sufficiently well to determine the steps that that 
approach implies for handling of the initial H in the 
model. But such understanding is demanded of most 
formulations of relativistic quantum theories for funda- 
mental reasons and the accuracy tests of similar kind 
can be consider a challenge for any scheme intended to 
solve the bound-state problem in theories with asymp- 
totic freedom. 



APPENDIX A: RGEP OF 4TH ORDER 

The RGEP weak-coupling expansion for terms of order 
g n , is written using terms of orders k < n, 



n—x 

G'n = E [/&,{(W)&-fc}' 



(Al) 



k=l 



If the interaction Hamiltonian Hi = H — Hq is propor- 
tional to g, the expansion in powers of g is the same as 
the expansion in powers of Hi . If Hi contains a polyno- 
mial in g with operator coefficients, then the expansion 
in powers of Hi is an intermediate step for obtaining 
an expansion for H\ in powers of g. Therefore, this ap- 
pendix provides a generic set of coefficients for expansion 
of £?/c(A) in powers of Hi. The argument A is omitted in 
what follows. The procedure of rewriting the expansion 
in powers of the bare g into the expansion in powers of 
gx requires a definition of g\ in the structure of H\ . This 
definition, in analogy to Eq. I|12|) , provides then a series 
expression for g\ (g) . This expression is inverted and sub- 
stituted into the formal expansion of H\ in powers of the 
bare g, producing the desired expansion in powers of g\. 

The first four terms in the perturbative expansion are 
written as 



[G 



[G 



m 1 n Hfji i H x r 



(A2) 
(A3) 

[=3] mn = bm\2nH m \H\iHin , (A4) 
H m lHi2H23H 3n . (A5) 



It is uderstood that the indices 1, 2, and 3, are summed 
over the entire range available for them in the theory. 

In order to write down expressions for the coefficients 
a, 6, and c, we introduce a set of auxiliary symbols. Their 
meaning becomes sucessively self-evident when one de- 
cifers them in the order they are given here. The vari- 
ables s and t have the meaning of the renormalization 
group parameter 1/A 2 . 



Fran F m E n , 

Gran — tPmn F n ^ 



(A6) 
(A7) 

Fran = fmn (E m ~ E n ) , (A8) 
G m ln = Gim + G\ n , (A9) 
F m ...lkn = F m ,,Jk + Fkn , (A10) 
A m l...kn\t) = fml—fkni (AH) 

Am...l+k...n (t) = A m ...i(t) A k ... n (t) , (A12) 



B m ...k(s) 



In this notation one obtains 



dtA, 



«(*)• 



(A13) 



(A14) 



b m 12n — [B m 12n+m2 — B m \2n] 

Em2 
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12 



E n 2 

G m ln G m i2 



F 



Bm2n-\-ml2 " J 



Gmln G n i2 
F n i2 



[B m 2n-\-nl2 

- B m2n ] (A15) 



ml 2 



6*2 



Cml23n 



E m l E m 2 

G 2 3n 
E n l E n2 

G132 
E\2 E m i 

G132 
El 2 E n i 
G\2m G 2 3 in 
F23m E m \ 
G\2n G23 n [ 
^23n E n i 
G\2m G13 n 
Fl3n E m i 
G\2n G\3 

in 

F\3m E n \ 
G\2m G\32 

F\32 E m \ 
G\2n G132 



[B m 321n — B m 321 n+lm + B nl 32i n+lm2\ 

\Bm\T3n ~ B m 123n+ln — -Bml23n+2n + B m l23n+ln2\ 
[B m 231n — -B r „231„+21 _ B m2 3i 1 + S m 231n+ml2] 
[B m 132n — B m 132n+21 ~ B m i32 n +nl + B m 132n+nl2] 

[Brn21n ^m21n+ml B m 21n-\-m32 + B m 21n+lm32\ 
B m \2n ~ B m 12n+nl ~ B m 12n+n32 + B m \2n+ln32\ 
\B m 2\n — B m 2l 31 + B m 21 n+ml3nj 

[B m i2n ~ Bml2n+nl BiYil2n-\-m31 ~\~ -^ml2n+m31nj 
[Bm21n — B 7n 21 n+ml B m 2]_n J rl32 H~~ -^m21n+m!32j 



Fl32 Enl 

G 



[B m 12n — B m i2n+nl — Brol2n+132 + -Bml2n+nl32] 



23/n 



E m 2 
Gniln G23 n 
E n 2 

G132 
-Ei 2 

Gmln G132 



Bmln Bni\n+m321 ^mln B m ln+m321 



+m321+m2 



-Fm321 
B m ln — B m \n+n32\ 
Fn32\ 



Frn32\ + -F2m 
Bmln — B m ln+n321+n2 



F n 321 + F 2r . 



B m ln ~ B m ln+m231 B m ln ~ B m ln+m2132 



F m 231 



F m 231 + F21 



E12 



Bmln B m ln+n231 B m ln B m .ln+n2 



132 



-Pn231 



-Fn231 + F21 



+ 



Gmln Gl2m G23m 

F23m 
Gmln Gl2n G23 n 

F23n 
Gmln Gi2m G132 

F132 
Gmln Gl2n G132 

F132 
Gmln Gl2n G13 

in 



Bmln Bmln+m21 B m ln B m ln-\-m321-\-m2 



F m 21 



F m 21 + F m 32 



Bmln Bmln J r n21 B m ln B m ln J rn321+n2 



F n 21 



F„21 + F n 



32 



B m ln ~ B m i n j rm 2i B m ln — B m ln+m2312 



F m 21 



F m 2i + -F231 



Bmln — B m \ n +n21 B rn ln — B m ln+n2312 



F n 21 



F„21 + F231 



Fi2n Ft, 



13m 



\Bmln B m ln+n21 Bmln-\-m31 ~t~ Bmln-\-m312n\ 



(A16) 



In relativistic applications, these formulae need a replace- 
ment of the differences of energies like E m — E n by the 
changes of invariant masses in the interaction vertices, 
and multiplication of G m n by the parent momenta P„ n 



in the vertices, see e.g. [TaLEfH]. 
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FIG. 6: The functions K{g x ) from Eqs. (Hi) to plot- 
ted for Wegner and RGEP cases, A = 2 GeV. In all cases a 
preferred value of g\ is clearly identified. These are given in 
Table I. 



TABLE I: The values of the couplings g\, which correspond 
to the minima of curves in Fig. [5] (the last digit is given with 
error margin of 5), and the corresponding results (obtained 
with the coupling constants given in this table) for the bound 
state energy \Ek\ (the exact value is 1), obtained by diago- 
nalization of W\ with A = 2 GeV, that was calculated in six 
successive orders of the weak-coupling expansion, k. 
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APPENDIX B: NUMERICAL DETAILS alytic formula from Appendix |X| 



Numerical integration of Eqs. © and © was per- 
formed using fourth-order Runge-Kutta procedure and 
an automated algorithm for generating expressions of or- 
der n from expressions of orders k < n. The numerically 
calculated matrix elements H\ mn were checked using an- 



The fits of the coupling constants g\ that were used in 
the perturbative evaluation of the windows W\ , were ob- 
tained using the functions K(g\) from Eqs. i|14|) and i|15|) 
that are plotted in Fig. The corresponding numbers 
for the coupling constants are given in Table I. 
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